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We investigate a possibility to restrict the chaotic inflationary scenario with large nonminimally 
coupled inflaton field <j> considered by Fakir and Unruh by means of the observed cosmic microwave 
background radiation(CMBR) anisotropy. This model is characterized by the condition £ > 1 and 
ip = 8ttG£4> 2 2> 1 where £ is the nonminimal coupling constant. We calculate the contributions 
of the long wavelength gravitational waves (GW) generated in the inflationary period to CMBR 
anisotropy quadrupole moment. We obtain the constraint A/£ 2 = 1.8 x lO -8 /^, where A is the self 
co uplin g and tpi means the initial value of ip. Combining this with previously obtained constraint 
y/ A/£ 2 w (<5T/T) rms = 1.1 x 10~ 5 , we conclude that the initial value has to be ip x fa 1.6 x 10 3 . If 
the self-coupling has a reasonably values of order 10 -2 , then £ « 10 4 and <f>i ~ 10~ 1 m p i where m v i 
is the Planck mass. 



I. INTRODUCTION 



Inflation is regarded as one of the most important conceptual progress in modern cosmology |l|,0]. This concept 
give us not only the solution of cosmological puzzles such as the horizon and flatness problem, but also the origin of 
the cosmological perturbations. Among various scenarios of inflation proposed so far, Linde's chaotic scenario seems 
to be a natural mechanism for the realization of the inflationary expansion [jj). In spite of many attractive features, 
the original scenario still have the fine tuning problem of the parameters. Namely one has to choose an unreasonably 
small self coupling constant, i.e. A < 1CU 12 , to have the correct magnitude of the density perturbation which is 
consistent with the observed tiny CMBR fluctuations. 

In order to overcome the difficulty, Fakir and Unruh proposed a new version of the chaotic model which has strong 
nonminimal coupling, i.e. £ > 1, to the scalar curvature |4|]. Note that we shall adopt in the present paper the sign 
convention for £ such that the conformal coupling means £ = —1/6. The chaotic inflationary scenario with nonminimal 
coupling has been investigated for various values of £ by Futamase and Maeda || . Their result does not exclude the 
strong coupling model of Fakir and Unruh. They have shown that the larger value of £ relaxes the fine-tuning problem 
of the self-coupling A, and this conclusion has been followed by Makino and Sasaki more rigorously |14|. 

Here we shall investigate Fakir-Unruh scenario from different point of view. Namely we shall study the effect of 
the scenario on the CMBR anisotropy. It is know that the long wavelength gravitational waves are generated during 
the period of the inflationary expansion, and they make a significant contribution to CMBR anisotropy through the 
Sachs- Wolfe effect [Q on lower multipoles [|ll],[R|. We will point out that Fakir-Unruh scenario has also a nonneg- 
ligible contribution to the spectrum of CMBR anisotropy by generated GW's, and then obtain the constraint on a 
certain combination of parameters by comparing the predicted quadrupole moment with the results of COBE-DMR 
observations. 

This paper is organized as follows. We first review Fakir-Unruh scenario and present the self-consistent inflationary 
solutions in Section 2. We then find the radiative modes of the metric perturbations in the spacetime generated 
by nonminimally coupled inflaton field, and derive the resulting spectrum of GW in section 3. Based on the result 
in section 3, we shall calculate the power spectrum of CMBR anisotropy and compare it with the observations of 
COBE-DMR in Section 4. Finally some discussions are made in section 5. We shall follow Misner, Thorne, and 
Wheeler [E2[ for the definition of Ricmann tensor, Ricci tensor and Ricci scalar. 
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II. BACKGROUND INFLATIONARY SOLUTIONS 



As explained in the introduction, we shall be in this paper interested in the chaotic inflationary scenario with large 
positive nonminimal coupling. Thus we shall take the following action as our model. 



S = 



x v-g 



R ^,2- 1 



2 k 



(2.1) 



where k 2 = 8ttG. Our definition of £ is the same as Fakir and Unruh, that is, c onformal coupling means £ 



Note that Futamase and Maeda used opposite sign for £. Varying the action (2.1), we obtain the field equations 



1 



Z<f> 2 Gfxu = (1 + 20M,u ~ 2£ 



1 



K" J \ 2 

- 2£0 (g^a - V^V*) + g» v V{4>), 
□</> + + V t<t> = 0, 



(2.2) 
(2.3) 



where □ = ~g 1 ' 2 d a (— g 1 ^ 2 g a ^dp) , and = dV/d(f>. Greek and Latin indices take 0, 1, 2, 3 and 1, 2, 3, respectively. 
For the spacetime we assume a homogeneous and spatially flat Robertson- Walker universe 



ds 2 = dt 2 -^{^dx'dx 3 . 
Equation( |2.2| ) yields the Hamiltonian constraint equation 



H 2 



3(1 + K 2 £(f> 2 ) 

where overdots denote time derivatives. The momentum constraint equation is 



H 



1 + K 2 £(f> 2 



The equation of motion for the inflaton field is obtained from (2.3) as 



3H<j) -6HH + 2H 



V,, = 0. 



Combining above equations, we obtain 



3H(f> 



k 2 £0 2 (1 + 6£) 



1 + k 2 £0 2 (1 + 6£) 

1 + k 2 £0 2 (1 + 6£) 1 
Using the "slow- roll approximations" such as 



4« 2 ^l/(0) - (1 + n 2 ^ 2 )V^\ 



« V (</>), 

the energy constraint and field equation take the following forms, respectively 



H 2 



3(1 + K 2 ^ 2 ) 

V(<P) 



2^ 



1 + k 2 ^ 2 (1 + 6£) 



{4 K 2 ^y(</»)-(i + K 2 e0 2 )^} 



1 



1 + k 2 £c/> 2 (1 + 6£) 



[4K 2 ^V(cf>) - (1 + n 2 ^ 2 )V^\ , 



(2.4) 
(2.5) 

(2.6) 
(2.7) 



(2.8) 

(2.9) 

(2-10) 
(2.H) 



(2.12) 
(2.13) 



These are the basic equations which determine the background solution. We shall consider the solutions of these 
equations in the case with the self-coupling constant and mass term, separately. 
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A. Chaotic Inflation by a self coupling 



Let us consider the chaotic inflation generated by the self-coupling constant of the inflaton. Thus we take 

v(4>) = \\^. 



Then the equations ( [2.12 ) and ( 2.13 ) give, respectively 



12(1 + K 2 £(/> 2 ) 

A0 3 



1 + 



8i 



l + K 2 £(/> 2 (l+6£) 



1 + k 2 £0 2 (1 + 6£) : 



(2.14) 
(2.15) 



Let us now define our basic variable ip = « 2 ?</> 2 which turns out to be very useful in our analysis. Since we consider 
the strong curvature coupling, let us consider the situation where 

i> > 1. 

This condition simplifies significantly the above equations as 

A-0 



i 



12k 2 £ 2 ' 



2A 



3k 2 £(1 + 6£) H' 



Therefore we have the following solution. 



H 



8£ 



1 + QC 
8£ 



■In 



l + 6£ 

Note that this solution is consistent with the slow-roll approximations. In fact we have 

< 1. 



i, 




8£ 1 


Hip 




1 + 6£V 



From the equations (2.16) and (2.19), we obtain a differential equation for the scale factor. 

1/2 



Mi 



12k 2 £ 2 



l 



8£ 



In 



(1 + 6$)^ \ai 



This equation (2.21) can be easily integrated as 

a(t) = a, exp [Hit - l{H t t) 2 ] . 

Thus we have the following Hubble parameter and scalar field. 

a{t) = a, exp [Hit - j(Hit) 2 ] , 
H(t) =Hi{\- 2 1 H i t) , 



where 



i>(t) = (i - 2-fHity , 



H? = 



7 



(2.16) 
(2.17) 

(2.18) 
(2.19) 

(2.20) 

(2.21) 

(2.22) 

(2.23) 
(2.24) 
(2.25) 

(2.26) 
(2.27) 
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These expressions indicate that the expansion rate is much more rapid than the rate of change in the Hubble parameter 
as well as the inflaton field. Since we are interested in the very long-wavelength tensor perturbations which give the 
largest contribution to the CMBR anisotropy today and such perturbations decouple at the early stage of inflation, 
we can approximate the above solution to the following simple forms: 



o(t) » a ie Bit , 
H(t) « Hi, 

Finally, we have to consider the amount of inflation sufficient to solve cosmological puzzles, 



a(t) 



exp 



exp 



ff *) =exp (/rH 

— §r #]=-g r (^-^/)>70 



(2.28) 
(2.29) 
(2.30) 



Fakir and Unruh commented ip f=s 1 signals the end of inflation 
gives 

V'i £ 80. 



(2.31) 

Thus the strong coupling, i.e. £ > 1, condition 

(2.32) 



B. Chaotic Inflation by a mass term 



Next, we consider a massive free inflaton field with the following potential 

V(cj>) = \m 2 4> 2 



Then the equations ( |2.12 )and ( 2.13 ) become, respectively 



H 2 



3H(p: 



6£(1 Ti^W) 
m 2 (j>(K 2 £,(f> 2 - 1) 



4£ 



l + s 2 ^(l + 6f) 



(i - 



1 + k 2 £0 2 (1 + 6£)' 
When ip ^> 1, these equations may be approximated as 

,2 



6£ V 1 + 6 ? 
2m 2 1 



(2.33) 
(2.34) 



V> 3(1 + 6C) H' 

We can find a scalar field grows exponentially so that the slow-roll approximations are not satisfied unless £ <C 1 



(2.35) 
(2.36) 







4£ 


Hip 




l + 2£ 



< 2. 



(2.37) 



Thus, we conclude the mass-term inflation does not occur in the model with strong coupling. This result is consistent 
with Futamase and Maeda @1. 



III. THE SPECTRUM OF GWS GENERATED DURING THE INFLATIONARY EXPANSION 



In this section we shall calculate the spectrum of GWs generated during the period od inflationary expansion in 
the framework of Fakir-Unruh scenario. 
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A. Radiative Solutions of Metric Perturbations 



First we derive the equation for the GWs as the linear metric perturbation. The method is exactly same with 
the case in general relativity (GR) |23|j . Namely we consider a small disturbance of the spatially flat Friedmann- 
Robertson- Walker metric g^ u : 

ds 2 = (g flu + h^ v )dx fl dx v . (3.1) 

The spatial flatness means that the metric perturbations can be expanded by a plane wave with a comoving wave 
number k : 

K(x, t) = [ d 3 kh x (k, Ty k - X et, (3.2) 

where e£f is the polarization tensor, and A = +, x is a mode of the polarization. We choose the synchronous gauge 
(/ioo = hoi = 0), and use the conformal time (a; = t), 

ds 2 = a 2 (r)dT 2 - {a 2 {T)5 l3 + hy) dx i dx j . (3.3) 



Then the linearized Einstein equations (2.2) lead to the radiative mode of the tensor perturbation: 

h k k = h\ k =81> = 0, (3.4) 

h'l(k, t) + (2- + h' x (k, t) + k 2 h x (k, t) - 0, (3.5) 



a 1 + %l) 

where dashes denote conformal time derivatives. The other components of the Einstein equations (|]^) are trivially 
satisfied with above conditions. 

Writing each component h\(k,r) of the GW perturbations as 

hx(k,r) = — !— Hx(k,T) = — =L==^ A (fc,T), (3.6) 
R{t) a-v/1 + ip 

we have the following equation for /x A . 

^(fc,r)+ ^fc 2 -^^ MA (fe,r)=0. (3.7) 

This is an ordinally "Schrodinger type" equation. If k 2 ^> \R" /R\, we can find a plane wave solution of h\ with its 
amplitude scaled by (ayT + ijij 

B. The Spectrum of GWs 

Now we will calculate the spectrum of the GWs derived above. We follow the previous works |Tl],[l2|,[l5| . Let us 
choose the normalization of a(r) as 

a(r) = — — - — vacuum re (— 00,— T2), (3-8) 

HiT 

= — X — radiation r 6 (t2, -^-), (3-9) 
r 2 

r 2 

— — matter t€(ti,to). (3.10) 
r o 

We assumed that the transitions across each phase are instantaneous and we matched a(r) and d(r) at the transition 
points. This will be an accurate approximation for the long-wave tensor perturbations we are interested in. Note that 
r is discontinuous across each transitions. 
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We shall calculate the amplitude of the GWs generated as qua ntum noises. The following consideration will help us 
to choose a convenient field variable to quantize. The equation (3/7) shows that each polarization state of the waves 
behave as a scalar field tp with a normalization factor of ^2k 2 s : 



2k 2 
1 + V 



(3.11) 



We omit a subscript A henceforth. This situation is analogous to GR case |L7j. We also neglected the other interactions 
of ip due to the mixture of ip' , which is a good approximation for our purposes. Thus we choose the scalar field ip to 
quantize as the same way in GR. 



(p(x,r) 



d s k 



a phasc (k)cp k (T)e tk - x + al hasc (k)cp* k (T)e 



ik X 



(3.12) 



where a p hasc and « Phase are the annihilation and creation operator at each phase(vacuum dominated, radiation dom- 
inated and matter dominated phases), and <^fc(r) is properly normalized solution of the equation (3.5) : 



¥>k(r) 



1 



2fc(27r) 3 / 2 a(r) 
1 



—iki 



—ikr 



2fc(27r) 3 / 2 a(r) 



fcr 



radiation. 



vacuum, matter, 



(3.13) 
(3.14) 



We used the fact that ip w const, in the inflationary phase (2.30), and ip — in the post inflationary phase. We also 
assumed our universe approaches the conformal vacuum mode function as t — > — oo, that is 



</?fc,VD 



-ikr 



2fc(27r) 3 / 2 a(r) 



for r 



-oo. 



(3.15) 



Where the quantities with "VD" means the quantities evaluated at the vacuum dominated phase, namely at the 
inflationary phase. The Bogoliubov coefficients relating operators at each phase are defined as : 



aRu(k) = d(k)ayr>(k) + c* 2 (k)a\j^{-k), 
&MD(fc) = c 3 (fc)avD(fe) + c| (—fc). 



(3.16) 
(3.17) 



where "RD" and "MD" mean "Radiation Dominated" and "Matter Dominated", respectively. Matching the field and 
its first derivative at the first transition: 



<p'k,vn(- T 2) = ci^ jRD (r 2 ) + c 2 tp*' KRD (T 2 ), 



we obtain 



ci 



C2 



"(fcr ) 2 



e 2ikT2 [1 - 2tkr 2 - 2(kr 2 ) 2 ] 



(3.18) 
(3.19) 



(3.20) 
(3.21) 



(fcro) 2 ' 

We can see these coefficients satisfy the Bogoliubov relations required for the orthonormality of mode functions <fk 



|ci| 2 -| C2 | 2 = l. 



(3.22) 



Since we are interested in GWs which are still well outside the horizon at the time of matter-radiation equality and 
will give the largest contribution to the CMBR anisotropy today, we may use simpler results up to second order 



ci 



c 2 



Hn 



(3.23) 
(3.24) 



G 



At the next transition the continuity conditions become 



Thus we obtain 



C3 



C4 



ci^,RD(y ) + c 2 ^,R D (y) = c 3 ^fc,MD(ri) + c 4 ^ )MD (ti), 

7~1 7~i 
Cl^fc,RD(y) + C2^fe !RD (y) = Cl^fc,MD( r l) + C2(p*fc !M D( T l), 



(3.25) 
(3.26) 



2(fcn) (fcr ) 2 
- 2ii 
H/k 



{1 - likT \ - 2(k Tl ) 2 } e lkTl/2 {1 - 2i/cT 2 - 2(fcr 2 ) 2 } e 2lA;T2 - e 3lfcri/2 



2(fcr 1 )(fcr ) i 



-3ifen/2 



{1 - 2i/cT 2 - 2(/ct 2 ) 2 } e 2lfcT2 - {1 + 2ikn - 2(/cti) 2 } 



2"l e -ifen/2 



with 



| C3 | 2 -| C4 | 2 = 1. 



Up to second order, we get simple solutions again as 



C3 



C4 



3iH 

2fc(fcT ) 2 ' 

2fc(fcr ) 2 ' 



(3.27) 

(3.28) 
(3.29) 

(3.30) 
(3.31) 



Since we assume each transition is instantaneous, the Universe will remain in the de Sitter vacuum. Thus we may be 
able to calculate the quantum mechanical two-point correlation function in the de Sitter vacuum. 



A (QM) = 



k 3 
(2ny 



d 3 xe ik x (0\<p(x,t)(p(0,t)\0) , 



where 



|0) = |de Sitter vacuum) 



(3.32) 



(3.33) 



Note that this assumption makes us to use the Heisenberg picture of quantum fields in which the states of vacuum 
don't evolve with time but the operators do Jis}| . For waves re-entering the horizon at the matter dominated era, we 
have 



a(q m ) = 



H 2 



2(2tt) 3 



3jl(fc7 

kr 



where ji{kr) is a spherical Bessel function of order one. According to our normalization, this is calculated as 



k 2 H 2 f3ji(A;r) 



kr 



where we used the following facts 



i v(A)c^(A')=«AV, 

ip{r > t 2 ) = 0. 



(3.34) 



(3.35) 



(3.36) 
(3.37) 



We match this result with our classical ensemble of GWs. The classical two-point function may be calculated by 
writing the classical GWs as 



h x (k,r) = A(k)xx(k) 



3ji(fcr) 
kr 



(3.38) 
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where x\ is a random variable with statistical expectation value 
Then the classical two-point function is 



A^\ xx ,=A 2 (k) 



1 2 



kr 



5\\> 



Comparing this with quantum two point function, we obtain 

k 2 H 2 8 



A\k) 



(2tt) 3 3tt 



where 



Since we know that H 2 = . n A t ' c2 , 



H 2 = — rrinl^v. 
3 p 



(16tt) 2 U 2 



(3.39) 



(3.40) 



(3.41) 



(3.42) 



(3.43) 



In the end, we find that the existence of the nonminimal coupling only affects the amplitude of spectrum via Hubble 
constant. 



IV. COMPARISON WITH OBSERVATIONS 

We shall now make a comparison between the theoretical prediction derived in the previous section with the 
observation of CMBR anisotropy. Using the Sachs- Wolfe relation, one can calculate the power spectrum of CMBR 
anisotropy generated by the long- wavelength GW fll| : 



<°*> = (E 



\ [7 + 2V /" 27r / ri 

a hn \ 2 \ = 367r 2 (2^ + l)j^jj J kdkA 2 (k)\F e (k)\ 2 , 



where the angle brackets denote averages over statistical ensemble of aj, and 

d ji (k(T - r)) 



Fi{k) 



dr 



d(k(r -r)) fc(r - r) 
je-2(kr) 2ji(kr) 



je+2(kr) 



_(2£- 1)(2£+ 1) (2£-l)(2£ + 3) (2£+l)(2£ + 3) 
We shall use the following numerical value for (a|) evaluated by White jl^] 

(a 2 ) = 7.74v. 

On the other hand, COBE-DMR group expresses this quantity in terms of Qrms-ps : 



(4) ^ ^ (1 



/rms-PS 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



Note that since this quantity Qrms-PS has already been handled st atist ica lly, we can ma ke a direct comparison between 
this expression with the theoretical prediction (a 2 ). Combining ( [3.43; ), (4.3) and (4A), we find 



\_ = (4tt)3 /80\ / Q rms - PS 
e 38.7 W'J V T 



(4.5) 
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According to COBE 4yr results |9|j2g], 

T = 2.728 ±0.004 K, (4.6) 

^5 rms — PS — 

18±1.4/LtK, (4.7) 
for n = 1 Harrison-Zel'dovich spectrum, thus we obtain the following constraint on the theory. 

*-«xio-(22), («) 



where we do not take the errors of the observations (4.6) and (4.7) into account, because our result (4.£) mainly 
depends on the initial value of fy. Furthermore, we have not considered the contribution from the scalar density 
perturbation which will be produced during early stage of the inflationary phase. Although a large fraction of the 
quadrupole anisotropy seems to be due to GW |lq] , we have to take that contribution into account for more precise 
discussions. Note that since the uncertainties just mentioned lea d to higher estimates of A/£ 2 , we may expect that 



the "real value" of A/£ is slightly lower than the obtained result (4.; 



Makino and Sasaki |14| derived the perturbation of spatial curvature lZ m on the hypersurfaces orthogonal to the 
matter rest frame : 



where N(tk) is a e-fold scale in which one is interested. lZ m is related to the gauge-invariant density perturbation as 



dp _ 2(w + 1) / k 



2 



., - , , M Km. (4-10) 
p ow + 5 \Ha J 

where w = p/p. Using our result ( |4.8| ), the density perturbation on the horizon scale (k = aH, N = 70) at a 
matter-dominated phase (w = 0) becomes 



5p 

P J HOR 




1.0 xlO- 4 ,/^. (4.11) 



The density perturbation is also related to the rms temperature fluctuation of CMBR via Sachs- Wolfe effect as 

=-h (4-12) 



T , 

where $ is the Newton potential, i.e. goo = 1 — 2$. On the other hand, $ is determined by the Poisson equation, 

k* = -A«G P ( 5 A=haH?( 5 A. (4.13) 



Since our interest is the density perturbation on present horizon scale, i.e. k = aH , we have 

(4.14) 



6T\ I /Sp 



HOR 




From the COBE-DMR result (6T/T) rms = 1.1 x 10 |20|, we obtain the initial condition of our universe as 

tpi fa 1.6 x 10 3 . (4.15) 

We must notice again, however, we have never paid attention to the ratio of tensor contribution to scalar one. Since 
smaller tensor contribution leads to smaller value of ipi by several factors, our discussion is reliable up to order of 
magnitude. 

Finally, let us consider the model of A = 10 -2 , which would be a reasonable value for some of particle physics 
model. Then the nonminimal coupling constant and the initial value of the inflaton field are constrained as 

£ « 10 4 , (4.16) 

A « 10"\ (4.17) 

rripl 

As Makino and Sasaki commented |Q , the relaxation of fine-tuning for A is nothing but a restatement of the issue of 
a large values of £. According to our results, furthermore, the initial inflaton field has to be smaller by one order of 
magnitude than that of chaotic scenario in framework of GR. Note that the smaller contribution of tensor perturbation 
leads to the smaller tpi and then smaller <pi/m p l. 







V. CONCLUSIONS 



We calculated the spectrum of GWs generated in the chaotic inflationary scenario with a nonminimally coupled 
inflaton under the condition where ip = K 2 ^(j) 2 ^> 1, and their contributions to the CMBR anisotropy via Sachs- 
Wolfe effect. Comparing the COBE-DMR 4yr observations with the predicted power spectrum of CMBR, we find 
the initial value ipi ~ 1.6 x 10 3 . Our result decreases a number of undetermined parameters in the theory, i.e. ipi, 
and indicates that if we choose A = 1CP 2 , £ w 10 4 and <pi/m p l ~ 10 _1 are required. In the chaotic inflationary 
scenario with minimally coupled inflaton, an argument based on the comparison between the energy content of the 
inflaton field and planck energy density places the upper bound of the initial value of the field. It would be hoped 
in the nonminimally coupled case that the same sort of argument will place an upper bound on the initial value for 
ipi. However it is the nonminimality that prevent us to have such a simple argument. In this context Hochbcrg and 
Kephart made some discussion on the energy density of a nonminimally coupled scalar field but their result 
seems to be not available in the present context because of the existence of self-coupling. It should be mentioned that 
the condition for Hubble constant to be smaller than the inverse of the planck time is easily satisfied in the above 
choice of parameters. In this sense Fakir and Unruh scenario is not excluded as a good candidate of the inflationary 
scenario. 

Finally we mention the chaotic inflationary scenario by induced gravity . The above analysis may be applied 
exactly same way and one has the same constraint of parameters when ip 3> 1. However this case has also another 
constraint on £ coming from solar system experiment on the omega parameter of Brans-Dicke theory since the theory 
is translated into Brans-Dick type. These two constraints are incompatible each other so that the chaotic scenario by 
induced gravity with the strong coupling seems not to work. 
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